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Introduction: There is ample evidence that Eu- 
ropa is currently geologically active [1-4]. Crater 
counts suggest that the surface is no more than 90 Myr 
old [5], and cryovolcanism may have played a role in 
resurfacing the satellite in recent geological times [6- 
10]. Europa’s surface exhibits many putative cryo vol- 
canic features [1,7-12], and previous investigations 
have suggested that a number of domes imaged by the 
Galileo spacecraft may be volcanic in origin [1,9] 
(Fig.l). Consequently, several Europa domes have 
been modeled as viscous effusions of cryolava [9,13]. 
However, previous models for the formation of silicic 
domes on the terrestrial planets contain fundamental 
shortcomings. Many of these shortcomings have been 
alleviated in our new modeling approach, which war- 
rants a re-assessment of the possibility of cryovolcanic 
domes on Europa. 



Figure 1. Some domes on Europa may have been 
emplaced as viscous effusions of cryolava onto the 
surface. The models presented here help constrain such 
cryovolcanic extrusions. 


New Modeling Approach: We have investigated 
the emplacement of cryovolcanic domes on Europa, 
exploring the effect of boundary conditions on the so- 
lution of the Boussinesq equation for pressure driven 
fluid flow in a cylindrical geometry. The continuity 
equation describing radial expansion of a Newtonian 
fluid with an unbounded (free) upper surface and a 
time-dependent viscosity is: 
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[14] found a similarity solution to the general form of 
(1) for a constant fluid volume with constant viscosity. 
This solution has been previously applied to the em- 
placement of putative volcanic domes on Venus and 
Europa [9, 15-16]. Here we offer an alternative simi- 
larity solution to (1) that eliminates the singularity at t 


= 0 inherent to the solution in [14] and allows for the 
investigation of associated plausible boundary condi- 
tions. This model also addresses the issue of time de- 
pendent changes in lava viscosity due to cooling. The 
similarity solution for flow thickness, h, is: 
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where ^ is a time transformation constant of the form: 
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The time constant that eliminates the singularity at t = 
0 is r = (3/4)^ (jilVf Vorolg- A variety of forms can be 
chosen for the time-dependent kinematic viscosity v(t). 
One plausible form is a cryolava viscosity that increas- 
es exponentially with time as it cools [17], v{t) = VoQ^'^. 
Fig. 2 shows the solution of a radially spreading, New- 
tonian fluid with Vo = 10^ m^/s (equivalent to 10^^ Pa-s 
for cryolava density of 10^ kg/m^) at four times. Here, 
the overall “shape” of the flow surface, as well as the 
aspect ratio at the final time, is very similar to the di- 
mensions of the domes in Fig. 1 [9] when r = 2.95 x 
10^ sec (34 days), and F =1 months. The total relaxa- 
tion time is ~ 3 years, which is well within the range of 
emplacement times for icy lavas predicted by conduc- 
tive and radiative cooling [9,18]. 
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Figure 2. Axially symmetric Newtonian fluid flow profiles 
obtained from (2). 


Skin Strength: Owing to water’s very low vapor 
pressure (-600 Pa), cryolava erupted into Europa’s 
low-pressure environment will boil violently until an 
approximately 0.5 m thick icy carapace forms; flows 
can then be maintained beneath this insulating layer 
[19]. For cryolava erupting at temperatures near 273 K, 
the latent heat of vaporization is approximately 7 times 
larger than the latent heat of fusion [19-20]. As a con- 






Figure 3. Growth rate of cryolava surface crust. Inset. A 
solid icy carapace 0.5 m thick will take ~ 3 days to grow. 

sequence, appreciable quantities of ice will be pro- 
duced for any given amount of lava that is boiled off, 
and conductive cooling will dominate [18, 20]. The 
growth of the icy crust atop the flow can therefore be 
modeled after [21]. Figure 3 shows an example of icy 
carapace growth as a function of time on a cooling 
cryolava at Europa’s surface conditions. 

If cryolavas on Europa are briny solutions [12, 22- 
23], their eutectic temperatures could be as low as 193 
K [24], and they may supercool 5-15° below their eu- 
tectics before crystallization commences [25]. The 
profiles in Fig. 4 show that after 3 years, approximate- 
ly 85% of a ~ 75 m thick flow would still be ductile, 
and hence able to relax and advance. The depth at 
which the glass transition temperature is reached in the 
flow can also be used as a gauge to determine the ex- 
tent of crustal growth [17, 26]. From Fig. 4 it is clear 
that the ~ 136 K glass transition temperature for water 
[27] is only reached down to shallow depths of ~ 3.6 
m within the flow. Hence the vast majority of the flow 
would still be mobile. 



laxed shape after 3 years (Fig. 2), only the top 10 m of the 
flow falls below 193 K, and only the top 3.6 m is below wa- 
ter’s glass transition temperature (i.e., 136 K). Hence, the 
flow will still be mobile. 

Nevertheless, crustal strength will have some ef- 
fect on bulk lava viscosity by providing a resistive 
pressure and slowing the flow of the fluid lava in the 
interior. The hydrostatic pressure of cryolava may ini- 
tially be large enough to entrain the icy crust, thereby 
ensuring continued dome relaxation [9, 19]. A similar 


scenario has been suggested for silicate lavas [28]. 
When taking into account the effect of the crust for 
terrestrial silicate lavas and cryolavas, the effective 
viscosity can be up to four orders of magnitude higher 
than the actual viscosity of the fluid [17-18, 29]. 
Therefore, the actual initial kinematic viscosity for the 
dome modeled in Fig. 2 may be significantly below 
10’ m’/s. 

If, based on their morphology, we assume that 
the domes in Fig. 1 are volcanic extrusions of cryolava, 
the projected kinematic viscosities for Europan lavas 
range from 10'^ - 10^ m^/s [9, 22]. The next step of this 
work will be to apply these models to other putative 
cryovolcanic domes imaged by Galileo to constrain the 
range of emplacement times and eruption rates for this 
particular style of cryovolcanism. 
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